bcftools variant calling